Complete break-up of the helium atom by proton and antiproton impact 
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We present a fully ab initio, non-perturbative, time-dependent approach to describe single and 
double ionization of helium by proton and antiproton impact. A flexible and accurate finite-element 
discrete-variable-representation is applied to discretize the problem on the radial grid in spherical 
coordinates. Good agreement with the most recent experimental data for absolute angle-integrated 
cross sections is obtained over a wide range of incident projectile energies between 3 keV and 
6 MeV. Furthermore, angle-differential cross sections for two-electron ejection are predicted for a 
proton impact energy of 6 MeV. Finally, the time evaluation of the ionization process is portrayed 
by displaying the electron density as a function of the projectile location. 

PACS numbers: 36.10.-k, 31.15.A-, 25.40.Ep, 25.43.+t 



^ : 

Oh, 
I 



Ph. 



> 
00 

o 

o\ 
o 



Ionization of helium by slow antiproton impact has re- 
ceived considerable attention in recent months, from exper- 
imentalists 111, and theorists alike 0, 0] • As a fundamen- 
tal, strongly correlated few-body collision system, it is at 
the heart of examining the reliability of the most advanced 
computational methods for atomic collision processes. Us- 
ing antiprotons has the major advantage of eliminating 
the complicated charge-exchange process that is competing 
with single ionization in the case of proton impact. 

Of particular interest regarding the validity of theoret- 
ical approaches is the low-energy region, where the pro- 
jectile speed \v\ is sufficiently slow that the Massey per- 
turbation parameter Q, where Zp is the charge 
of the projectile, is becoming larger than unity. Note 
that \Zp\/\v\ = 1 for incident (anti)protons with energy 
Ei ~ 25 keV. Consequently, perturbative treatments based 
on one or even a few terms of the Born series are entirely 
inappropriate for the problem at this and lower energies. 
Instead, non-perturbative treatments, either based upon 
the close-coupling expansion with an appropriate set of 
basis functions 0, 0] or the direct solution of the time- 
dependent Schrodinger equation (TDSE) on a numerical 
space-time grid 0, |3i S| are required. The latter method 
is often referred to as the time-dependent close-coupling 
(TDCC) approach. 

The interaction between the charged projectile and 
the target can also be viewed as a temporal ultrafast 
half-cycle-like pulse, which is responsible for the (multi- 
ple) ionization process of the target The full width 
at half maximum of such a pulse can be estimated as 
2^/2>b/\vp\. For (anti)proton impact, this is of the order 
10 attoseconds for "intermediate" energies around 100 keV. 
Furthermore, the peak strength of the electric field Wmax 
is of the order Zp/\h\. Even for (anti)protons as the light- 
est "heavy" particles, Z^max is approximately two atomic 
units at an impact parameter |b| ~ Ioq, where ao = 
0.529 X 10~^° m is the Bohr radius. This results in peak 
intensities of the half-cycle pulse as high as lO^^W/cm^. 
Consequently, charged-ion impact presents an interesting 
alternative to intense laser-pulse techniques in the study of 
ionization processes in strong electromagnetic fields within 
an ultrashort time window 91. 



Various experimental datasets for single and double ion- 
ization of helium by antiproton impact were published over 
the years. While there is generally good agreement at 
high incident energies, both between different experimen- 
tal datasets and predictions from various theoretical mod- 
els, the situation is much less clear for incident energies of 
20keV and below. Most recently, Knudscn et al. [l| pub- 
lished results for single ionization of helium by antiproton 
impact that differed substantially from those obtained ear- 
lier by Andersen et al. [l^l and Hvclplund et al. lT| . While 
the recent experimental data of Knudscn et al. l| were re- 
produced fairly well by various non-perturbative calcula- 
tions 0,0,0] I discrepancies of 15% or more still remained. 

Here we report results from a fully ab initio numeri- 
cal study of the helium break-up problem by antiproton 
and proton impact over a wide range of incident energies 
between 3 kcV to 6 MeV, as well as the corresponding to- 
tal and double-differential cross sections (DDCS). Our ap- 
proach, whose validity is not restricted to a particular pro- 
jectile energy range, thus provides a unique opportunity 
for studying the multiple ionization dynamics induced by 
charged ions from weak to strong perturbations. The con- 
figuration space of the target electrons is discretized via a 
finite-element discrete- variable representation (FE-DVR). 
This highly flexible and accurate grid-based approach com- 
bines the numerical advantages of basis-function expan- 
sions in small intervals with an easily adaptable spatial 
grid to account for the radial dependence of the electronic 
wavefunction close to and far away from the nucleus. 

The dynamical response of the system is obtained by 
propagating the initial wavepacket, defined on the DVR 
gridpoints, through a recently developed time-dependent 
Arnoldi-Lanczos algorithm [ij, ■ The collision system 
is governed by the time-dependent Hamiltonian 
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The essential complexity, compared to previous work in- 
volving a spatially homogeneous laser field [l^ , is the rep- 
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resentation of the two-body Coulomb interaction 



2Ai + l [ri,i?(t)]; 



between the projectile and the target. Here ri and r2 
are the coordinates of the two helium electrons, R{t) 
is the coordinate of the projectile, and [a;,y]<(>) = 
min(max){x, ?/}. All these coordinates are defined relative 
to the He2+ ion, which is fixed at the origin. 

We use a straight-line trajectory R{t) = b + (dp — vt) 
for the incident projectile starting at do at an impact pa- 
rameter b. This is the only physical approximation made 
in our treatment. It should be sufficiently accurate even 



at a projectile energy as low as 3 keV, which is the low- 
est energy considered in this work. At this energy and 
|5| ~ 1 ao, which is the impact-parameter regime with the 
largest contribution to the cross sections, the scattering an- 
gle in the laboratory system is estimated to be merely 1°. 
This results in a small relative momentum transfer | g | / 1 | 
of 1.8%, where Pi is the projectile's initial momentum. 

As a result of using a classical trajectory for the pro- 
jectile, the total angular momentum i, magnetic quan- 
tum number M, and parity 11 of the collision system are 
no longer conserved quantities that could be taken advan- 
tage of in a fully quantal partial-wave expansion. However, 
there is still one conserved quantity, namely the reflection 
symmetry of the electronic wavefunction with respect to 
the collision plane. The latter is defined by the trajectory 
of the incident projectile and the impact parameter b. 

We explicitly build this reflection symmetry into our FE- 
DVR wavefunction by writing 
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Instead of ordinary coupled spherical harmonics 
ytih {''^i'''^^), we introduced the angular basis 
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These basis functions are normalized according to 
{QlX\QvX') = 5i^v,5i,v^5ll'5mm'- Using the reflection 
symmetry relates the wavefunctions for the magnetic quan- 
tum numbers M and —A/, in the same way as it does 
for electron-impact excitation and coplanar ionization and 
ionization-excitation processes [IJ] . Consequently, we only 
need to include M > in Eq. ([3]) and thus the size of the 
problem is reduced to nearly half of what it would be with- 
out adopting the above symmetry. Taking advantage of the 
reflection symmetry, we use 55 angular partial waves gen- 
erated by setting (/i, -^^)max = (3,3,3,3) to expand 
the wave function in Eq. This is equivalent to the 101 
partial waves used in Ref. ^31- 

To calculate the angle- and energy-integrated cross sec- 
tions, 399 DVR grid points were set up in a spatial box of 
''max ~ 80 Op , while a smaller step size and 799 points were 
employed for the DDCS. We also extended the box size to 
160 ao in this case to ensure converged results. 

The total cross section (angle- and energy-integrated) is 
obtained as 
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where P{b,Ei) = J \{^kik2\^{t))\'^dkidk2 represents the 
probability for single or double ionization at fixed values 
of |6| and Ei. Figure [T] exhibits our results for single ioniza- 
tion of helium by antiproton impact for projectile energies 
between 3 keV and 5 MeV. In the theoretically most diffi- 
cult low-energy regime, we obtain excellent agreement with 
the most recent experimental data of Knudsen et al. [ij. 
For projectile energies of 20 keV and above, our results are 
slightly lower than those of Foster et al. Q and thus in 
better agreement with experiment. On the other hand, 
most experimental data near the cross section maximum 
around 100 keV and beyond lie above both our results and 
those of Foster et al. Having performed extensive conver- 
gence checks, we are confident in the numerical accuracy 
of our predictions and currently have no explanation for 
the remaining discrepancies. 

Figure [2] depicts the corresponding results for the double 
ionization process. Although the size of the experimental 
error bars [2] is substantial and thus limits the meaning of 
comparing the absolute numbers between theory and ex- 
periment, we note that our results are in excellent agree- 
ment with those of Foster et al. [3| in the limited range of 
projectile energies where their data are available. 

To gain further insight into the dependence of the joint 
two-electron response on the sign of the projectile charge, 
we show in the left panel of Fig. [3] an example for 100 keV 
antiproton impact on helium at an impact parameter of 
0.5 ao, while the right panel shows the corresponding re- 
sults for proton impact. We first note that not much hap- 
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FIG. 1; (Color online) Cross section for single ionization of 
helium by antiproton impact. Experimental data obtained at 
CERN by Andersen et al. [T^] (CERN90), Hvelplund et al. [lit 
(CERN94), and Knudsen et al. [l[ (CERN08) are compared 
with TDCC 01 and the present FE-DVR predictions. 
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FIG. 3: (Color online) Radial electron density after antiproton 
(left panels) and proton (right panels) impact on helium at an 
energy of 100 keV for an impact parameter of 0.5 ao. Starting 
from an initial distance of —50 ao, the positions of the projectile 
shown in the snapshots, from top to bottom, are +10, +20, and 
+40 ao relative to the center of the target, respectively. 



FIG. 2: (Color online) Cross section for double ionization of he- 
lium by antiproton impact. The experimental data obtained at 
CERN by Andersen et al. [3 (CERN90), Hvelplund et al. [ll| 
(CERN94), and Knudsen et al. 0] (CERN09) are compared 
with TDCC and the present FE-DVR predictions. 



pens until the projectile is already 10 oq passed the target 
(top panels). For this case, and also a little later when the 
projectile is 20 ao behind the target (center panels), the ra- 
dial electron densities show a substantial double-ionization 
component, whose signature is a significant density at large 
values of both ri and r2. The most energetic electrons have 
moved to about 25 and 50 oq, respectively. Finally, when 
the projectile is 40 oq beyond the target (bottom panels), 
we see the characteristics of both double and single (charac- 
terized by significant densities when only ri or r2 is large) 
ionization processes developing. By this stage, the most 
energetic electrons have moved as far away as 100 oq from 
the He^+ center. When comparing the electron densities 
for the two projectile charges, we recall that the antipro- 
ton cannot capture any of the electrons and in fact pushes 
them away, whereas the capture channel may be important 
for proton impact. While the single- ionization signals for 
the two projectiles resemble each other, the major differ- 



ence between the results for the two projectiles occurs in 
the double-ionization region of ri « ^2. Apparently the 
proton is trying to attract at least one electron and hence 
causes a reduction in the probability for both electrons 
moving out with the same speed. 

In order to portray the two-electron emission in a more 
quantitative way, it is important to consider angle-resolved 
cross sections, e.g., the DDCS for two-electron ejection 
without observing the electrons' individual energies. This 
particular DDCS is obtained as 



— ^ 27r / bdb / dki / dfca 
dkidk2 Jo Jo Jo 



2 



(6) 



where cr; denotes the Coulomb phase and {ki, ^2) is 

the partial-wave amplitude in momentum space. 

Figure |4] exhibits such DDCS results for proton im- 
pact double ionization of helium for an incident energy of 
6MeV. The results are for the coplanar geometry, where 
the momentum transfer and the momentum vectors of the 
two ejected electrons all lie in the same plane. We com- 
pare our FE-DVR predictions for cuts with fixed values of 
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blc ionization. Not surprisingly for such a high projectile 
energy, their ratio of ~ 400 is consistent with what one 
would expect from the shake-off mechanism. 

In summary, we have investigated the complete break-up 
problem of helium by proton and antiproton impact by us- 
ing a time-dependent non-perturbative FE-DVR approach. 
This is a prime example of a highly correlated four-body 
Coulomb process, whose description remains a major theo- 
retical and computational challenge. Generally good agree- 
ment with the latest sets of experimental data for both in- 
tegrated and differential cross sections was obtained. At 
lower projectile energies (< 20keV), our antiproton results 
clearly show that the cross section for double ionization 
decreases with decreasing projectile energy. For the angle- 
resolved DDCS presented here, converged results require 
to account for contributions from ejected electrons with 
energies of several hundred eV, whereas the angular de- 
pendence is essentially established by low-energy electrons 
with energies of less than ~ 25 cV. 

We thank Prof. Hclgc Knudscn for sending us experi- 
mental data in numerical form and helpful discussions and 
Dr. James Colgan for clarifying problems with Ref. [1]. 
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grid allocation TG-PHY090031. 



FIG. 4; (Color online) DDCS for proton impact double ioniza- 
tion of helium for an incident energy of 6 MeV, a fixed detection 
angle 9i for one of the electrons, and a variable detection an- 
gle 02 for the second electron. The experimental measurements 
of Refs. [Tsl . [l6| were normalized to the converged DDCS at the 
large peak for 6i = 5°. 

9i that were generated by Foster et aim from the original 
experimental data of Fischer et al. [l5| after the analysis 
by Schulz et al. 0. 

Note that only ejected electrons with energies up to 
25 eV were observed experimentally, while converged re- 
sults for the DDCS defined in Eq. ^ require energies up 
to ~ 400 eV to be counted in this case. Although the mea- 
sured signal thus only corresponds to about one third of 
the converged DDCS, we see that the shape of the DDCS 
is essentially determined even with the cutoff at 25 eV. 

We see good agreement with the experimental data and 
clearly reproduce the systematic shift of the large peak 
to the right with increasing values of 9i. Furthermore, 
the angle between 9i and the maximum in the DDCS as a 
function of 02 decreases slowly from about 120° for 6i — 5° 
to about 100° for di = 55°. Similar results were published 
by Foster et al. ^ , but both the graphical presentation and 
the theoretical magnitudes given in Fig. 2 of their paper 
contain errors 17|. 

Finally, we obtained values of 3.54 and 8.97 x 10"^ Mb, 
respectively, for the total cross sections for single and dou- 
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